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Abstract 

The influence of the heating mechanism on the kinetic energy densities of the components of 
a vibrated granular mixture is investigated. Collisions of the particles with the vibrating wall 
are inelastic and characterized by two coefficients of normal restitution, one for each of the two 
species. By means of molecular dynamics simulations, it is shown that the non-equipartition of 
kinetic energy is not affected by the differential mechanism of energy injection, aside the usual 
boundary layer around the wall. The macroscopic state of the mixture in the bulk is defined by 
intensive variables that do not include the partial granular temperatures of the components. 

PACS numbers: 45.70.-n,47.70.Nd,51.10.+y 
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A fundamental issue for non-equilibrium physics is the identification of the variables 
needed to characterize a given macroscopic state [l]. Granular gases are a widely studied 
example of out of equilibrium systems ^. Due to the inelasticity of collisions, they do 
not have an equilibrium state and, in order to keep the granular medium fluidized, external 
energy must be continuously supplied in some way. An extension of the usual hydrodynamic 
equations for molecular gases hasproven to provide an accurate description of granular 
gases in many different contexts |3|. This description involves the granular temperature, 
defined from the mean kinetic energy of the particles, just by analogy with the definition of 
temperature in kinetic theory of molecular fluids. 

For granular mixtures, one of the consequences of inelasticity is that the energy equipar- 
tition required in equilibrium systems is not verified. The granular temperatures of the 
components of a mixture, defined from the average kinetic energy of each species, are differ- 
ent. This feature was pointed out many years ago {4], and a systematic study of it has started 
in the last years. Once established that the partial temperatures of the several components 
of the mixture are different, some important conceptual questions arise: should they be 
incorporated into the set of intensive variables needed to identify a given macroscopic state 
of a granular mixture, replacing the single temperature field employed in the description of 
one-component systems? If the answer is negative, are the partial temperatures determined 
by the usual set of hydrodynamic fields (including the temperature of the mixture) and 
the properties of the two types of particles? It is important to realize that for ordinary 
mixtures, the partial temperatures of the components are not included in the hydrodynamic 
description. Moreover, kinetic theory then implies that the temperature fields of all the 
components are the same, i.e. equipartition applies locally 

In the theoretical study of granular mixtures, it is usually assumed that the temperature 
fields of each of the especies are not needed to describe the system j6|. In a kinetic theory or 
statistical mechanics description, this implies that the temperatures of the components can 
be determined from the other macroscopic parameters, as indicated above, since the (normal) 
distribution function of the system is assumed to be determined by the hydrodynamic fields, 
so all the properties of the system turn out to be functionals of them 

Recently, Wang and Menon [7] reported the results of some event-driven simulations of a 
granular mixture and reached the conclusion that the heating mechanism affects nonequipar- 
tition of energy, even in the bulk of the system. This would imply that the details of the 
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driving used to inject energy into the system can not be ignored in describing mixtures of 
inelastic gases even far away from the energy source. If this were actually the case, the 
generalized hydrodynamic-like equations for granular gases should include the partial tem- 
peratures of the components, and the existing derivations of them should be deeply revised. 
In this paper, the influence of the differential boundary heating on the violation of equiparti- 
tion is reconsidered by means of molecular dynamics simulations. It is found that the details 
of the driving mechanism do not affect the relations between intensive quantities aside from 
a boundary layer next to the energizing wall. 

In the event-driven simulations whose results will be described below, the system is 
composed of two kinds of smooth inelastic hard disks, labeled 1 and 2. Both types of disks 
have the same diameter a, but different masses, mi and m2- The number of particles of each 
species is A^^i and A''2, respectively, being the total number of particles = A^^i + A'"2 = 420, in 
all cases. There is an external gravitational field acting on the system, so that each particle 
i is subjected to a force —rriigoez, where go is a positive constant and is the unit vector 
in the direction of the positive z axis. The inelasticity of collisions is modeled by constant 
coefficients of normal restitution. There are three of them: an, 022, and ai2 = 021, where 
aij refers to collisions between a particle of species i and a particle of species j. They are 
defined in the interval < a^- < 1. For the simulations reported here, ai2 has been chosen 
as given by ai2 = (an + a;22)/2. The system is open on the top and its width is W = 70a. 
The latter value has been chosen small enough so that the steady state with gradients only 
in the vertical (z) direction is stable, avoiding the development of transversal instabilities 

n n n 

[8|, 19|, llOj. To eliminate undesired boundary effects induced by the lateral walls, periodic 
boundary conditions in the horizontal (x) direction have been employed. 

To keep the system fiuidized, energy is continuously supplied through the wall located at 
the bottom {z = 0) of the system. The specific way in which this wall is modeled has been 
chosen from a compromise between several considerations. The physical issue of interest, 
unequally heating of the two species, should be isolated, avoiding additional effects on the 
properties of interest induced, for instance, by the oscillations of the wall. Also, the heating 
should be, as much as possible, the (idealized) model of some experimental process. In this 
context, "thermal" walls in which the velocity of a reflected particle is not related with its 
incident velocity, are not at all neutral for the "adjustment" of the dispersion of the velocities 
of the two species, since they destroy all the velocity correlations by definition. In addition. 
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it does not seem clear to which type of experimental boundary condition for granular fluids 
they correspond to, if any. Thus a vibrating wall with a sawtooth velocity profile is chosen, 
so that all the particles colliding with the wall find it with the same upwards velocity vw 
Moreover, the wall is supposed to move with an amplitude much smaller than the mean 
free path of the particles in its proximity. Consequently, an accurate description is obtained 



by considering it as fixed at 2 = 
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12l |. Therefore, the dynamics of this wall does not 
induce any additional space or time dependence. The mechanism for the differential heating 
must be now introduced. It will be attached to the difference in the coefficient of normal 
restitution, awi and aw2, for the collisions of the two types of particles with the vibrating 
wall. When a particle of species i having a velocity v collides with the wall, its component 
Vx remains unchanged, while the z component is instantaneously modified to 

v'^ = - {1 + awi){v^ - Vw)- (1) 

Note that a collision of this kind is only possible if < 0, since the wall is treated as located 
at 2; = 0. From Eq. ([T]), it is obtained that 

-vl = {1 + awi) {vw + \vz\) [(1 + o:wi)vw - (1 - avy)!^^^!] ■ (2) 

Therefore, when a particle collides with the vibrating wall, its kinetic energy decreases as a 
consequence of the collision if 

1 + awi /„x 

\Vz\>Vc=- Vw (3) 

1 - awi 

When the average velocity of the particles of one of the species approaching the vibrating 
wall is larger than Vc, this species loses energy in the collisions with the wall, on the average. 
Of course, the net balance of energy flux through the vibrating wall when considering both 
species must be positive since, in the steady state, energy must be continuously supplied by 
the wall to compensate for the dissipation in collisions. Note that this effect, cooling of one 
of the species by the vibrating wall, is not possible when the particles collide with the wall 
in an elastic way, as Vc diverges in this case. On the other hand, it can occur when thermal 
walls are considered, as made for instance in [3]. 

Four differential drivings will be considered in the following: (1) awi = 0.999 and aw2 = 
0.5, (2) awi = 0.995 and aw2 = 0.9, (3) awi = 0.9 and aw2 = 0.995, and (4) awi = 0.5 
and aw2 = 0.999. In Fig. [1] the density, ni{z), and temperature, Ti{z), profiles for each 



of the two species in the steady state are shown for each of the above drivings. In all the 
simulations reported in the figure, an = 022 = 0.93, m2 = 5mi, and Ni = N2. These are 
the same values as considered in ref. Q]. The velocity of the vibrating wall, vw, has been 
chosen in each case such that the system remains in the dilute regime. The specific values 
used for each of the above drivings are: vw / y/gocr = 7 for driving (1), 4 for driving (2), 5 
for driving (3), and 6 for driving (4). To illustrate how the mechanism of energy injection 
actually affects to each of the two species, the average variation of the kinetic energy of the 
particles when colliding with the vibrating wall in the steady situation, < Acj >, has been 
computed for the two kind of particles. The value of the ratio mi < Ae2 > /m2 < Aei > 
for each of the four drivings specified above is: 0.19 for driving (1), 0.47 for driving (2), 0.85 
for driving (3), and —6.20 for driving (4). The disparity in the above values, clearly shows 
that the input of energy per unit of mass for each of the components depends very strongly 
on the particular differential driving being used. The negative sign for driving (4) is due to 
the partial cooling of species 1 upon colliding with the vibrating wall, that corresponds to 
the possible scenario pointed out above. 

In Fig. 1, it is seen that the profiles of the hydrodynamic fields are different for the 
different drivings, everywhere in the system. Actually, the same happens if the density 
ratio, n2{z)/ni{z), or the temperature ratio, 7(2) = T2{z)/Ti{z), are plotted. This is not at 
all surprising, since the hydrodynamic profiles and also the density and temperature ratios 
are expected to depend on the boundary conditions. In monodisperse systems with an 
elastic vibrating wall, it has been shown that v^r plays the role of an scaling factor for the 
amplitude of the hydrodynamic fields of the gas in the steady state 13(]. Mathematically, 
this property follows by realizing that the energy injected by the wall into the system is 
$VK = NmgQVw/W, and it is a consequence of the inelastic Navier-Stokes hydrodynamic 
description. From this point of view, the verification of the scaling in a region can be 
considered gnature of the bulk of the system. 

For the two-component system with differential heating we are considering, it can be ex- 
pected on physical grounds that the relevant quantity characterizing the vibrating boundary 
at a macroscopic level is again the injected energy flux, and that in the low density limit 
it scales the hydrodynamic fields in the same way as vw does in the one-component case. 
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FIG. 1: (color online) Number density, rii, measured in units of u"^, and granular temperature, Tj, 
measured in units of rriiago, profiles for each of the the two species of particles, along the vertical 
direction. Height is measured in units of a. The masses of the two kinds of particles are related 
by 771-2 = 5mi. There are 210 particles of each species and the width of the system is W = 70a. 
Four levels of differential heating at the vibrating wall, labeled from 1 to 4, are shown; they are 
specified in the main text. The density profiles for species 2 (red online) exhibit higher maxima 
than those of species 1 (black online). Each temperature profile for species 2 (red online) is always 
above the profile for specie 1 (black online) corresponding to the same driving. 

Then, dimensionless temperature and density fields are defined by 

< = ^, (4) 



and 



■w 

respectively. Consistently, an scaled coordinate z* is introduced by 



(6) 



The scaled density and temperature profiles in the new length scale for the same systems 
as in Fig. [1] are plotted in Fig. O The measured values of for each of the four drivings, 
measured in units of mgQ^'^a~^^'^, are: 74.62 for driving (1), 63.67 for driving (2), 82.83 for 
driving (3), and 81.00 for driving (4). Although the curves do not overlap perfectly, it seems 
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FIG. 2: (color online) Dimensionless number density, n*, and granular temperature, T* profiles, 
defined in Eqs. ^ and ([S]), for the two species of particles along the vertical direction, in the 
dimensionless length scale introduced by Eq. ([6]). The system and the several differential drivings 
are the same as in Fig. [TJ 

fair to conclude that there exist well defined scaled profiles for each of the components, aside 
from a boundary region next to the vibrating wall. Actually, the observed discrepancies can 
be due to finite density effects, similar to those found in p^], since the maximum value of 
total number density reaches values n = rii + n2 ^ 0.1, which are not very low. In this 
context, it is worth to insist that the scaling for one- dimensional systems only holds in the 
dilute limit. 

Now that the appropriate scaling to compare the hydrodynamic profiles has been identi- 
fied, it can be used to check the hydrodynamic nature of 7, in the sense that, in the bulk of 
the system, it is determined by local properties of the hydrodynamic fields of the mixture. 
In Fig. [31 the profile of 7 is plotted as a function of the dimensionless scale z*. It is observed 
that the curves corresponding to the four differential heating mechanisms collapse in a wide 
region of the system, actually in most of it since the fluctuations for large values of z* are 
due to the small number of particles present there. This collapse is a very strong indication 
of the hydrodynamic character of 7, and it shows that the departure from equipartition is 
determined by the local properties of the mixture, being unnecessary to introduce indepen- 
dent temperature flelds to describe the behavior of each of the two species. It is important 
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FIG. 3: (color online) Temperature ratio profile, 7, for the two species of particles along the 
vertical direction. The system and the differential heatings are the same as in Fig. [TJ The vertical 
z* coordinated is measured in the dimensionless units defined in Eq. ([B]), which involve both the 
parameters defining the system and the particular heating being used. 

to realize the relevance of the scaling when investigating the bulk properties of 7. Although 
trivially it does not affect the ratio, the value of z* following from a given z depends, not 
only on the parameters defining the system itself, but also on the differential heating being 
used, since both determine the amount of power injected through the vibrating wall in the 
steady state. It can be wondered why the collapse of the 7 profiles is clearer than the col- 
lapse of the partial temperatures showed in Fig. [2J A probable reason is that finite density 
effects are more relevant for each of the component temperatures than for its ratio. The 
same kind of collapse was obtained for other systems with mass ratio between 0.5 and 5, and 
with values of the coefficients of normal restitution up to roughly 0.8. For smaller values of 
these coefficients, the one-dimensional hydrodynamical profiles become unstable and a more 
involved analysis is needed. 

In summary, the simulation results reported here show that when energy is supplied to 
a granular mixture through a vibrating wall, the details of the collision mechanism for each 
of the species do not affect the extent of non-equipartition, beyond the usual boundary 
layer. This means that a macroscopic description does not require to introduce the partial 
(granular) temperatures of the components. We expect a similar conclusion to hold for 
other kind of differential boundaries and also for sheared granular fluids. Of course, a 
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different question is the theoretical prediction of the departure from equipartition, measured 
for instanc e by the partial temperature ratio. The theoretical studies carried out up to 
0, [is], 16] refer to ho mog eneous systems, with the exception of the dilute tracer limit 



now 



in a vibrated system 



13, 



18| . although the results have been compared to experimental 



data obtained in spatially inhomogeneous vibro-fiuidized systems, and a relative qualitative 
agreement has been found 
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